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Abstract 

We analyze and demonstrate, numerically and experimentally, the self-healing effect in 
scaled propagation invariant beams, subject to opaque obstructions. We introduce the signal 
to noise intensity ratio, a semi-analytical figure of merit, explicitly dependent on the features 
of the beams and the obstructions applied to them. The effect is quantitatively evaluated 
employing the Root Mean Square deviation and the similarity function. 
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1. Introduction 


An optical beam that is subject to a partial obstruction propagates showing certain degree of 
self-reconstruction in the obstruction domain. This phenomenon is usually referred to as self- 
healing (SH) of the beam. The SH effect has been demonstrated and studied mainly in propa¬ 
gation invariant beams (PIBs), as Airy f 11, Bessel (2j|5), Caustic |[6j, Mathieu and Weber Q 
beams. 

There are other interesting beams whose transverse intensity profile is also invariant under 
propagation, changing only their scale. We refer to these beams as scaled propagation invariant 
beams (SPIBs), to distinguish them from the PIBs. Examples of SPIBs are the Hermite-Gauss 
(HG), Ince-Gauss (IG), and Laguerre-Gauss (LG) beams, which are solutions of the parax¬ 
ial scalar wave equation in three different coordinate systems fl4|p3| . Considering that these 
beams show a kind of propagation invariance, it is also interesting to investigate their self¬ 
reconstruction capabilities. In this paper we analyze, demonstrate and evaluate the SH effect in 
SPIBs. 

In section 2 we first recall the structure of SPIBs. Then we propose a method to analyze 
theoretically the SH effect in a generic SPIB, introducing the signal to noise intensity ratio, a 
semi-analytical figure of merit of the self-healing process, explicitly dependent on the features 
of the beam and the obstruction applied to it. This figure of merit, which is computed in base of 
analytical diffraction formulas, provides an indirect theoretical measurement of the degree of 
self-healing that can be expected for an obstructed beam. In addition to this indirect theoretical 
figure of merit, there is necessity of computing the true degree of self-reconstruction of prop¬ 
agated beams. This quantitative assessment is performed by means of the Root Mean Square 
(RMS) deviation and the similarity function (S) 0. In section 3 we demonstrate and evaluate 
by means of numerical simulations and experimentally the SH in HG and IG beams. 


2. Theory 

2.7. Structure of scaled propagation invariant beams 

The complex amplitude of the SPIBs (HG, IG and LG beams) includes a common factor with 
the modulation of a Gaussian beam, which is characterized by the Rayleigh length zo- Other 
parameters, expressed in terms of zo, are wo = (Azo/^) 1//2 , the beam radius waist, w{z) = 
wo[l + (, z/zo ) 2 ] 1 ^ 2 , the beam radius as a function of z, and R(z) = z[l + (zo/z) 2 ], the curvature 
radius of the quadratic phase. The Gouy phase, which is proportional to arctan(z/zo), will be 
specified for each one of the beams. Next we describe the structure of SPIBs. Lor brevity we 
only present the analytical expressions for the HG and the IG beams. 

The HG beam is the solution of the paraxial Helmholtz equation expressed in rectangular 
coordinates. Its complex amplitude is given by © 


HGB(x,y,z) (^*) H n (f^y) exp 



— ikz + i(f) (z) — ik 


2 R(z) 


( 1 ) 


where Hj{j = /, n) is the j-th order Hermite polynomial, r = \/x 2 + y 2 is the radial coordinate, 
and (j)(z) = (/ + n+ 1) arctan(z/zo) is the Gouy phase. The complex transmittance of the HG 
beam, at the waist plane (z = 0) is expressed as 

HGB(x,y,z = 0) = E 0 H, (&x) H n (^y) exp (-^) . (2) 

The other SPIB to be considered, the IG beam, is the solution of the paraxial Helmholtz 
equation expressed in elliptical cylindrical coordinates. Such coordinates (£,rj,z) are defined 
by the relations v = /(z)cosh(<^)cos(rj), y = /(z)sinh(^)sin(rj), and z = z, where f(z) = 







fow(z)/wo is the semifocal separation dependent on z and fo is the semifocal separation at the 
waist plane. From these formulas it is established that the domains of t, and T] are [0,°°) and 
[0,27 t), respectively (19]. The complex amplitudes of the even ( e) and odd (o) IG beams are 
given by 


IGBlJt; ,77 ,z)=cjfcCZ(i?, e)C™ ( 77 , e) exp 
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IGB;^,r 1 ,z) = s^ S™{it, , e)Sp(ri ,e)exp 
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(4) 


where C™ and S™ are the Ince polynomials 
(m, p) that are in the ranges 0 < m < p and 


w(z) 

of order p and degree m. The integer indices 
m< p for even and odd functions, respectively, 


must have the same parity, i. e. (—1 ) p ~ m = 1. The elliptic ity p arameter of IG beams is £ = 
2/q/wq, and the Gouy phase is (j)(z) = (p + 1) arctan(z/zo) 151. At the waist of the IG beams, 
their complex amplitudes reduce to 


IGB e p 


,m( ( 5) T 7>^ = 0) = c C™ (it;, £)C™ (ri , e) exp (—^) , 


(5) 


IGB^ m (!;,7 ] ,z = 0)=sS™(it;,e)S™(T 1 ,e)ex V (-^), 
In Eqs. {Tg, Eq, c and ^ are normalization factors. 


( 6 ) 


2.2. Analysis of self-healing in scaled propagation invariant beams 

To analyze the SH effect in SPIBs let us consider the setup depicted in Fig. □ n is assumed that 
a SPIB of complex amplitude b(x,y), which arrives from the left side of the setup to the plane 
z = 0, is partially obstructed at this plane. 



Fig. 1. Optical setup: a beam of complex amplitude b(x,y) partially obstructed at the plane 
z = 0 is transformed into b Q (x,y). Then, b 0 (x,y) propagates to the distance z becoming 
b p 0 (x,y). 

The complex amplitude of the obstructed field at the plane z = 0 can be expressed as 

b Q (x,y ) = b(x,y)[\ - o{x,y)\, (7) 

where o(x,y ) is a binary function, equal to 1 at the obstruction area, and 0 otherwise. Let us 
denote as b p (x,y) the field that would propagate to the distance z if the beam b(x,y ) were not 
obstructed. This field can be expressed as 


b p {x,y)=FR z {b{x,y)} 


( 8 ) 























where FR Z is the operator that represents Fresnel free propagation to a distance z. The field that 
propagates to the distance z when the obstruction is present, can be expressed as 

b%{x,y)=FR z {b 0 (x,y)} = b p (x,y)-n(x,y), (9) 


where 

n(x,y) = FR z {b(x,y)o(x,y)}. (10) 

According to Eq. the non-obstructed propagated field b p (x,y ) appears as a part of the ob¬ 
structed propagated field. The other function in Eq. ([9]), n(x,y ), plays the role of a perturbation 
of the non-obstructed field. To analyze the SH effect it will be helpful to know the structure of 
n(x,y). This can be done applying in Eq. (10 ) the formula for the Fresnel transform of a product 
of functions (D, to obtain 


«(*>?) = {E*(x,y,z)b p (x,y)} ® O , 


(ID 


where 0(u,v) is the Fourier transform of the obstruction pupil o(x,y), <S> denotes the convo¬ 
lution operation, and E{x,y) = exp[(/;r/Az)(.* 2 -Fy 2 )]. According to Eq. ([9]), a good approxi¬ 
mation of the non-obstructed field b p (x,y) will be recovered in the propagated field b%(x,y) if 
the intensity of n(x,y) is negligible in comparison to the intensity of b p (x,y). Next we establish 
estimated values of such intensities and propose a criterion for self-reconstruction and the prop¬ 
agation distance required to fulfill it. To determine the intensities we require optical powers of 
the involved fields [in Eq. ([9])] and the areas where these fields are distributed. We first establish 
the transverse widths of the field b p (x,y) as 


Ax b = 2a x w(z ), 

A y b = 2cCyw(z ), (12) 


where a x and 0Cy are constant parameters and w(z) is the radius of the beam Gaussian factor, at 
the propagation distance z. To establish the widths in Eq. and the corresponding parame¬ 
ters a x and CCy , we specify the beam limits (along any horizontal or vertical line) at the positions 
where the edge intensities reduces by a factor of Xje 1 (respect to the peak intensity). The hor¬ 
izontal field width is the maximum of the widths obtained for the different horizontal lines. A 
similar criterion provides the vertical field width. It is also important to establish the transverse 
widths Ax 0 and A y 09 of the function 0(x/Az,y/Az) that appears in Eq. (11). As particular case 
we assume that o(x,y) is a square obstruction, of width a , for which we obtain 


Ax 0 = Ay 0 = 2Xzja. (13) 

Considering general features of the convolution operation in Eq. fTT] ) we estimate that the 
transverse widths of the perturbation field n(x,y ), along the horizontal and vertical axes, are 
respectively 


Ax n = Ax b + Ax 0 , 

Ay n =Ay b +Ay 0 . (14) 

The areas of the domains of b p (x,y) and n(x,y), in terms of the widths defined in Eqs. © 
and (D are given by 


A b = Ax b Ay b , 
A n = Ax n Ay n . 


( 15 ) 





On the other hand, we denote the optical powers of b p (x,y ) and n(x,y ) as P b and P n . Considering 
that such powers are invariant during propagation, we can compute them at the plane z = 0. 
The powers P b and P n results from integrating \b(x,y)\ 2 in the whole plane z = 0 and in the 
obstruction domain respectively. Therefore, we can establish the average intensities I b = P b / A b 
and I n = P n /A n , which correspond to the field b p (x,y) and the perturbation function n(x,y ), 
respectively. Now we can introduce the parameter 


■L n 


(16) 


a positive quantity, which is referred to as signal to noise intensity ratio. An important result, 
proved below, is that the ratio g is limited by an upper bound, explicitly dependent on the 
parameters of the beam and the obstruction applied to it. Introducing the relative power P nb = 
PnjPb , Eq. (16) can be expressed as 

A y=QPnb- (17) 

Ab 

Considering Eqs. © to ( |15| ) and performing some algebra, Eq. ( |17| ) leads to the relation 
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= \ QaPnb + at-a, 
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(18) 


where a = 0Cy/a x , a p = 
obtain the inequality 


(a + l)/2 and aj = (a — l)/2. Considering that Ax 0 /Ax b > 0, we 


a l- a l 

OtP n b 


< 2 , 


(19) 


that provides a minimum bound for Q. On the other hand, considering definitions for Ax 0 and 
Ax b , Eq. ( p~8| ) can be transformed into 
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( 20 ) 


which corresponds to the propagation distance required to obtain a given signal to noise inten¬ 
sity ratio Q. Now, since the expression under the square brackets in Eq. ( [20] ) must be positive, 
one obtains the inequality 


Q< 


1 

CCP n b 



7lW 0 

aa x 


2 
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( 21 ) 


that gives a maximum bound for g, which complements the minimum bound, given by ©.it 
is noted that the maximum value of Q, in < [2T| >, is dependent on the obstruction size a, the relative 
power P nb , and the beam features related to the parameters a, a x , a p and As a consequence 
of this result, given the parameters that specify the beam and the obstruction applied to it, there 
is a limit for the parameter g, which corresponds indirectly to a limit for the degree of SH that 
can be attained on propagation. When the beam is not obstructed both parameters a and P nb are 
null and the upper bound for g becomes °o. In any other case, the upper bound value of g is 
finite. 












3. Numerical and experimental results 

To evaluate the validity of the quantitative formulae proposed in the previous section we study 
the SH of HG and IG SPIBs. For the experimental generation of the beams under test we em¬ 
ploy amplitude computer-generated holograms (CGHs), displayed in a twisted nematic liquid 
crystal (TNLC) spatial light modulator (SLM) G3- The holograms are designed to include 
the obstruction as a feature of the generated fields, simplifying the experimental set-up. Let 
us briefly describe the features of these CGHs. We assume that the CGH is used to encode 
the optical field s(x,y) = a(x,y)exp(z' 0 (%,y)) where the amplitud a(x,y) is a normalized posi¬ 
tive function and the phase <j>(x,y ) takes values in the range [—7i,7i\. The transmittance of the 
amplitude CGH that allows the generation of the field s(x,y) is given by 

h(x,y) = hb(x,y)+a(x,y) cos [ 0 (*,y) - i2n(uox + voy)] , ( 22 ) 

where (wo, vo) are the spatial frequencies of a linear phase carrier and h b (x,y) is a background 
function that makes positive definite the function h(x,y) (17) . In order to propitiate a high signal 
to noise ratio in the generation of the field s(x,y) we must chose a background h b (x,y) with low 
power and low bandwidth. In the CGHs that we implement, h b (x,y) is chosen as a soft Gaussian 
function. 

The Fourier spectrum of the CGH in Eq. ([22]) is 

H(u, v) = H b (u, v) + ^S(u + wo, v + v 0 ) + ^S*(-u - w 0 .-v- v 0 ) (23) 

where H b (u,v) is the Fourier transform of h b (x,y), S is the Fourier transform of s(x,y), and 
S* is its complex conjugate. We assume that S(w,v) corresponds to one of the desired SPIBs 
(HG or IG beams). To isolate the SPIB, a band-pass filter, centered at frequency coordinates 
(—wo, —vo), is placed at the Fourier domain of the CGH. 

The experimental setup designed to generate the SPIBs, is depicted in Fig. [2] In this setup, 
an expanded and collimated He-Ne laser beam is used to illuminate the CGHs codified on 
a TNLC-SLM. Linear polarizers P\ and Pi, orthogonal to each other, are required to obtain 
mostly amplitude modulation in the SLM. The Fourier transform of the amplitude CGHs is 
generated by the lens L 3 . The field S(u + wo, v + vo), consisting in one of the SPIBs (HG or 
IG beams) appears at the open pupil in the spatial filter (SF) plane. For convenience, a dark 
area corresponding to the required obstruction, is optionally encoded in the SPIB domain. The 
intensities of the generated fields are recorded with a CCD camera at different distances along 
the propagation axis. 


He-Ne Laser ]| 


Ll 

Fig. 2. Experimental setup to demostrate self-healing of HG and IG beams. These fields, 
that include optionally a dark area (or obstruction area), appear at the first order in the 
Fourier domain of the CGH encoded in the TNLC-SLM. The CGH Fourier transform is 
generated by the lens L 3 . The fields are recorded with a CCD camera at different distances 
z. 


Pi SLM P 



In Fig. [3] we show the intensities generated numerically of a HG beam (a) and an odd IG 
beam (b) that will be employed to illustrate the SH effect in SPIBs. The parameters for the HG 







beam are / = n = 8 , and for the odd IG beam are p = 8 and m = 2. The intensities of the beams, 
experimentally generated with amplitude CGHs, using the optical setup in Fig. [2] are shown 
in Fig. 3(c,d). The focal length of the lens L 3 was / = 75cm, and the waists of the generated 
beams, at the plane of the SF, were 200/1 and 137/i respectively. 



Fig. 3. (a, b) Numerical and (c, d) experimental intensities of HG and IG beams, employed 
to illustrate the SH effect in SPIBs. 


In order to test the SH effect in partially obstructed SPIBs, the CGHs used to generate the 
fields in Fig. [3j were modified to include different dark square areas (obstruction domains), at 
the center of each field. The intensities of the obtained fields, at different propagation distances, 
are displayed in figures |4] to [ 6 ] For each figure the results corresponding to the numerical sim¬ 
ulations and the experiments appear at the top and the bottom images respectively. For each 
horizontal array of images, in any of the figures, the field at z =0 (at the left extreme of the 
array) presents the obstructed domain in focus. The widths of the fields of view, in the different 
images, are proportional to the width w(z) of the Gaussian factor, at the corresponding distance 
z. 


In Fig. [4] we present the results for the HG beam, with indices l = n = 8. The waist radius 
and obstacle width in this case are respectively w 0 = 200 jum and a = 2.25 w 0 . Employing the 
definitions in section 2.2 we obtained the parameters a x = OCy = 3 and P n j, ^0.1. For these 
parameters we employed Eqs. ( [T9] ) and ( [21] ) to obtain the allowed values of Q in the range 
[9.90:18.72). For each Q in this allowed range we can compute the required distance z/zo, us¬ 
ing Eq. ^2A\ . For example, if we chose the ratio Q = 16.5, the computed normalized propagation 
distance obtained is z/z^ = 1.2. It is noticed that the upper value of Q in the interval is avoided, 
since in this case z diverges to infinity. Both, the numerical and experimental results in Fig. 0 
correspond to normalized propagation distances (z/z 0 ) in the range [0 : 1.2]. The propagation 
distance for the fields in Fig.[4]marked as (d) and (h) correspond to the signal to noise intensity 
ratio <2=16.5. The results obtained for the same HG beam employing an increased obstruction 
area of size a = 4.16 w 0 , are displayed in Fig. [5] In this case, the parameters a x and OCy remain 
unchanged and the relative powers take the increased value P n b = 0.25, obtaining the allowed 
values of Q in the range (4.03 : 6.32). The different recorded fields (numerical and experimen¬ 
tal) correspond to propagation distances z/z 0 in the range [0 : 1.6], with incremental step 0.4. 
In particular the propagation distance z/z 0 = 1.6 corresponds to the signal to noise intensity 
ratio Q = 5.94. This Q value, and the associated propagation distance, correspond to the fields 
displayed in the images marked as (e) and (j). 

In the case of the IG beam, the SH demonstration is performed with an obstruction of width 
a = 2.7 w 0 (where w 0 = 137 gm). Others parameters of the IG beam are the constants a x = 2, 
OCy = 3, and P n t = 0.358, for which we obtained the allowed values of Q in the range (2.79 : 
7.02). The results (numerical end experimental) displayed in Fig. [6] correspond to propagation 
distances z/z 0 in the range [0 : 1.28] with approximated steps of 0.43. The propagation distance 
z/z 0 = 1 -28, in this range, correspond to the intensity ratio Q = 5.34 (smaller than the limit value 
Q = 6.16). In the numerical and experimental results, displayed in figures [4] to [6] it is noticed 








Fig. 4. Numerical (top images) and experimental (bottom images) results for self-healing 
of a HG beam subject to a medium sized square obstruction, at the propagation distance 
z/z 0 : (a,e) 0, (b,f) 0.4, (c,g) 0.8, and (d,h) 1.2. 



Fig. 5. Numerical (top images) and experimental (bottom images) results for self-healing 
of a HG beam subject to a relative large square obstruction, at the propagation distance 
z/z 0 : (a,f) 0, (b,g) 0.4, (c,h) 0.8, (d,i) 1.2, and (e,j) 1.6. 


a relatively good self-reconstruction of the field inside the obstruction domain. In contrast, the 
field outside the obstruction area has been clearly affected by the perturbation term n(x,y). This 
qualitative observation is verified by the quantitative assessment of the propagated fields, in the 
next subsection. 


3.1. Quantitative evaluation of self-healing 


In section 2.2 we introduced the parameter Q that represents a rough form of signal to noise 
ratio in the self-reconstruction of the desired field h p , during propagation of the obstructed field 
bo. As a complement of this semi-analytical assessment of the self-healing process, here we 
evaluate this process employing two figures of merit: the Root Mean Square (RMS) deviation 
and the Similarity (S) function |j4j. 



























(e) (f) (g) (h) 


Fig. 6. Numerical (top images) and experimental (bottom images) results for self-healing 
of a IG beam subject to a small sized square obstruction, at the propagation distance z/zfr 
(a,e) 0, (b,f) 0.43, (c,g) 0.85, and (d,h) 1.28. 


The RMS deviation of the intensity of the obstructed field, I p = \b p Q \ 2 , respect to intensity of 
the non-obstructed field, I p = \b p \ 2 , is given by: 


RMS = 


J (P - pl§) 2 dA 

W 

WYdA 


(24) 


where dA is the differential area, w is the domain where the RMS is evaluated, and /3 is a 
constant that allows the best fitting of the intensities I p and I p , obtained from the relation, 
dRSM/dp = 0. On the other hand, the similarity is defined as: 


c_ m,m) 

m\ my 


(25) 


where (\b%\,\b p \) = f w (b p ) (b p )*dA is the inner product of two fields, meanwhile \b p \\ = 
L fw(bo) ( bo)*} l ^ 2 dA and \ \b p \\ = [J w (b p ) ( b p Y] l / 2 dA are the norm of the fields, and aster¬ 
isk (*) denotes complex conjugation. For each propagation distance (z), the RMS and S are 
computed in a square domain Q of size a(z) = a{w(z)/w 0 ), referred to as field internal domain. 
The scale change in this domain is equal to that of the beam itself, during its free propagation. 
Additionally, we evaluated numerically the RMS and S for a domain external to Q . 

Values for the RMS deviation and the similarity S in the domain Q, for the HG beams in Fig. 
3(a)| subject to different centered square obstructions, and propagated to different distances, are 
displayed in Fig.[7Ja,b). In addition to the medium and large obstruction, which correspond to 
the ones employed in the numerical and experimental results (for the HG beam) just discussed, 
we also considered here a small obstruction, of width a = 1.68 w 0 . Similar results for the IG 
beam in Fig. |3(b)| employing square obstructions of two different widths, are displayed in 
Fig.[7jc,d). In this case the medium and large obstructions have widths 2.80 w 0 and 4.16 w 0 , 
respectively (with w 0 = 200/lm). As noted in the results of Fig. [7] the behavior of the similarity 
is opposite to that of the RMS deviation. It is also noted that the RMS metric presents more 
sensible variations along the propagation ranges. An interesting fact is that the RMS attain its 
best (minimum) value at certain propagation distances. E.g. the RMS for the HG beam with 











the medium sized obstruction shows a minimum at z/z 0 ~ 2.3. Another case of the IG beam, 
with the medium sized obstruction that shows a minimum RMS at z/z 0 ~ 2.2. This behavior is 
consistent with the prediction in section |Z2| of an upper bound value for the signal to intensity 
ratio Q. 



(a) 



(c) 



(b) 



(d) 


Fig. 7. Behavior of RMS deviation (a,c) and Similarity function (b,d) within the internal 
domain Q, at different propagation distances of HG (a,b) and IG (c,d) beams subject to 
different obstruction areas. 


To complement the results in Fig.[7J we also computed the RMS and S for the already consid¬ 
ered beams and obstructions employing the domain external to Q (within the area of b p (x,y)). 
It is interesting to note that the results, displayed in Fig. [8j show a clear degradation of the 
external propagated fields. Fortunately, this degradation also saturates at certain propagation 
distances. 


4. Final remarks and conclusions 

We have analyzed the SH effect in SPIBs, a kind of propagation invariant beams with scale 
change. In particular we demonstrated numerically and experimentally the effect in the cases 
of HG and IG beams. 

We established a semi-analytic parameter to asses the degree of SH, corresponding to the 
signal to noise intensity ratio ( Q ). This quantity is computed in terms of parameters that repre¬ 
sent the features of both the beam and the obstruction applied to it. We obtained that for a given 
set of these parameters, there exist a restricted range of values for Q. Moreover, we presented 
a formula to compute the required propagation distance for each Q value in the allowed range. 
This analysis establishes indirectly that the degree of SH in SPIBs can not be achieved beyond 
certain limit, regardless the propagation distance of the obstructed beam. Such restriction is 
confirmed qualitatively in the images of self-reconstructed beams, obtained numerically and 
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Fig. 8. Behavior of RMS deviation (a,c) and Similarity function (b,d) outside the domain 
£1 at different propagation distances of HG (a,b) and IG (c,d) beams subject to different 
obstruction areas. 


experimentally, in section section 3 (figures [4] to [6]) . 

To complement the semi-analytic assessment of SH by means of the parameter g, we eval¬ 
uated numerically the effect by computing the RMS deviation and the Similarity for the ob¬ 
structed propagated beams, respect to the non-disturbed beams. We first computed such pa¬ 
rameters in a domain Q (referred as internal domain), which corresponds to a scaled version 
of the applied obstruction. The RMS deviations computed in this domain, present a reduction 
during propagation, up to a limit that is dependent on the beam and obstruction features. The 
appearance of this limit roughly provides a quantitative confirmation of the limit in the SH, that 
was analytically established by means of the g ratio. Although the similarity S is less sensible 
that the RMS, the SH limit can also be noted at close views of the S plots (which for brevity are 
not shown). 

We also computed the RMS and S in a domain external to Q (within the area of the beam 
b p (x,y)). The results (in Fig. [8]) are interesting, although not surprising. They show that the 
RMS increases and the similarity reduces at the the external domains. 

The analytical discussion regarding the signal to noise intensity ratio (g) is an original con¬ 
tribution in this paper that can be extended to analyze the SH in other propagation invariant 
beams (Bessel, Airy, etc). Another contribution has been the separate evaluation of the RMS 
and similarity in an internal and an external domain. The results in this context improves the 
understanding of the SH effect. 















































